Universality in chaotic quantum transport: 
The concordance between random matrix and semiclassical theories 
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Electronic transport through chaotic quantum dots exhibits universal, system independent, prop- 
erties, consistent with random matrix theory. The quantum transport can also be rooted, via the 
semiclassical approximation, in sums over the classical scattering trajectories. Correlations between 
such trajectories can be organized diagrammatically and have been shown to yield universal answers 
for some observables. Here, we develop the general combinatorial treatment of the semiclassical di- 
agrams, through a connection to factorizations of permutations. We show agreement between the 
semiclassical and random matrix approaches to the moments of the transmission eigenvalues. The 
result is valid for all moments to all orders of the expansion in inverse channel number for all three 
main symmetry classes (with and without time reversal symmetry and spin-orbit interaction) and 
extends to nonlinear statistics. This finally explains the applicability of random matrix theory to 
chaotic quantum transport in terms of the underlying dynamics as well as providing semiclassical 
access to the probability density of the transmission eigenvalues. 
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Closed mesoscopic systems with sizes between the 
atomic and macroscopic possess statistically different en- 
ergy spectra depending on whether the corresponding 
classical system is regular or chaotic P, 0) ■ A semiclassi- 
cal approach to such systems, valid in the effective limit 
of ?i 0, leads to trace formulae where the density of 
energy states is approximated by sums over the classical 
periodic orbits of the system which form stable fam- 
ilies for regular systems while being unstable and isolated 
in chaotic ones. A hallmark of the energy statistics is 
the form factor, a two-point correlation function approx- 
imated by a double sum over periodic orbits. By pairing 
orbits with themselves for chaotic systems or members 
of their families for regular ones, the difference between 
their corresponding energy spectra can be directly linked 
to the properties of the underlying dynamics 0, |g| . 

For (quantum) chaotic systems, there is the further 
conjecture [3] that the statistics of the energy spectra are 
universal (depending just on the symmetry of the system) 
and identical to those of the eigenvalues of large random 
matrices 8] , originally employed to model the spectra of 
atomic nuclei. However, the semiclassical pairing of pe- 
riodic orbits with themselves ^ only led to agreement 
with the leading order term of the random matrix the- 
ory (RMT) form factor. Recently, additional correlated 
periodic orbits were identified, treated and shown to pro- 
vide exact agreement with RMT for short times [ol. flo|. 
This involves orbits which come close to themselves in 
an 'encounter', whose occurrence can be estimated using 
the ergodicity of the classical motion, and partner orbits 
which can be constructed, due to the local hyperbolic- 
ity, to cross the encounter differently. For long times the 
correlations remain unknown, but the form factor can be 
obtained through resummation of short orbits 11- 13|. 




FIG. 1. (Color online) (a) The trajectories for M2 form a 
closed cycle if the dashed trajectories, which contribute with 
negative action, are traversed backwards, (b) To contribute in 
the semiclassical limit, the trajectories must be nearly iden- 
tical apart from in small encounter regions which they can 
traverse differently, (c) By untwisting the encounter, the di- 
agram in (b) can be redrawn as the boundary walk of a tree. 



systems, obtained by attaching scattering leads, as seen 
for example in an experimental study of the electronic 
transport through quantum dots [l3l- Theoretically, we 
start with the transmission subblock t of the scattering 
matrix connecting asymptotic states in the (two) leads. 
The transmission eigenvalues of the matrix T = tH, and 
their moments M„ = Tr[r"] relate to the electronic flow 
through the system. For example in the low temperature 
limit the first moment is proportional to the conductance 
MM- 

For ballistic chaotic systems, modeling the scattering 
matrix by a random matrix from the circular ensembles 
was proposed and shown to be consistent with a diago- 
nal semiclassical approach 17| . For the low moments Mi 
and M2, all off-diagonal contributions were evaluated in 



201 ■ while the calculation of general M„, but only 
for the first several off-diagonal terms, were performed 
221 . In all cases, the results agree with RMT. 
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Hallmarks of the underlying dynamics persist in open 



The purpose of this letter is to exhibit the mathematical 
reasons behind this agreement and establish the general 
equivalence between semiclassics for open systems and 
RMT of the circular ensembles. The derivation extends 
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FIG. 2. (Color online) The graph in (a) and its boundary walk represent a trajectory quartet which would contribute to M2. 
Following steps 1-5 we first cut the top encounter node (keeping 02 connected to 12) and arrive at (b). As oi is now attached 
to a node of degree 6 we insert a new link to obtain (c). Reconnecting the link from 12 to 02 we obtain a second trajectory 
quartet in (d) which exactly cancels the contribution from (a). Performing steps 1-5 on (d) reverses the chain to recover (a). 



to all three main symmetry classes: unitary for chaotic 
systems without time reversal symmetry (TRS), orthog- 
onal for systems with TRS and symplectic for systems 
with spin-orbit interaction. It also remains valid for the 
nonlinear moments. 

Our approach holds in the universal regime where the 
dwell time, the average time spent inside the system, 
is much longer than the Ehrenfest time te, the time 
needed for a wavepacket of size Ap, the Fermi wave- 
length, to grow to the system size L and delocalize. Un- 
der chaotic dynamics, with Lypunov exponent A, we have 
Te ~ A^^ In (L/Ap) and when no longer small compared 
to the dwell time, RMT stops being applicable. However, 
such Ehrenfest time effects have been incorporated into 
the semiclassical framework for all diagrams at leading 
order and some subleading order diagrams for low mo- 
ments [2^-[2^. Our systematic approach may then be 
useful beyond the universal regime. 

RMT results. RMT provides the joint probability dis- 
tribution of the transmission eigenvalues [27| which can 
be integrated to obtain the transport moments, as was 
performed for the conductance and its variance 28, 2^. 
Other quantities were limited to diagrammatic expan- 
sions [30[ until the connection to the Selberg integral was 
explored 3l|, |34| . Since then there has been much inter- 
est and success in calculating the moments Af„ from the 
circular ensembles [33l - l36l |. 

Semiclassical diagrams. Semiclassically, the elements 
toi of the scattering matrix are approximated l^, 37, 3^ 
by a sum over the trajectories 7 which start in channel 
i in one lead and end in channel o of the other. They 
contribute a phase exp(iS'.^/?i) with their action so 
that Mn is approximated by a sum over 2n trajectories 
of which half contribute with positive action and travel 
from channels ij to Oj while the other half contribute 
with negative action and travel from channels ij+i to Oj 
(we identify in+i with ii). Geometrically, if we reverse 
the direction of the trajectories with negative action, the 



trajectories would form a single cycle visiting zi, oi, Z2, . . . 
in turn, as in Fig. [TJa). The phase involving the actions 
oscillates in the semiclassical limit unless the total ac- 
tion difference is small on the scale of h. To obtain the 
statistical properties of the moments we average over a 
range of energies so that oscillating phases wash out and 
only trajectory sets which achieve this small action dif- 
ference contribute consistently. These, as for closed sys- 
tems, come close in encounters while being nearly iden- 
tical elsewhere (in 'links'), as in Fig.[ljb). 

For Ml (conductance), we have trajectory pairs start- 
ing and ending together, and the contributing semiclassi- 
cal diagrams were identified and treated [l9| precisely by 
cutting open the periodic orbit diagrams of the form fac- 
tor. For M2, which is related to the shot noise, the con- 
tributing trajectory qua drup lets can be formed by cut- 
ting periodic orbits twice [20(] , and this approach has been 
generalised in the complementary work of Ref. [39| . 

The contribution of each semiclassical diagram is a 
simple product of its constituent parts |20l| . If the leads 
carry A^i and iV2 channels respectively (with a total 

= iVi -I- A'2) then each incoming channel provides the 
factor A^i and each outgoing channel N2. More impor- 
tantly, each link provides a factor oi 1/N and each en- 
counter —N. Assuming A^i ^ iV2, the order in N^^ of 
each diagram is the difference between the number of 
links and the number of encounters (and channels). The 
leading order diagrams can then be redrawn as trees, 
or rather as paths around the tree along the so-called 
'boundary walk' 21|. The encounters are untwisted to 
become roundabout nodes, the links edges and the chan- 
nels leaves, so Fig. [IJb) morphs to Fig.[ljc). The leading 
order of all Mn was obtained in ^21] by recursively gen- 
erating the trees. Higher order diagrams involve closed 
cycles and a graphical representation provided moment 
generating functions at the next two orders [22| which 
match an asymptotic expansion of RMT results [36[ . We 
now build on this to show exact concordance between 
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FIG. 3. (Color online) The trajectory sextet in (a) form a single cycle or the permutation (123). Untying the node with 
leaves 03 and 02 (keeping 03 connected to 13) breaks the single cycle into two and corresponds to multiplying (123) by (23) 
to obtain (12) (3) as in (b). Repeating the untying we move through (c) and (d) until we reach separated hnks in (e) or 
the identity permutation. Inverting the steps we can represent the original diagram in (a) as the primitive factorization 
(12)(13)(13)(23)=(123). 



semiclassics and RMT for all moments and to all orders. 

Cancellations. First we show that the vast majority 
of possible semiclassical diagrams cancel. Since each en- 
counter leads to a minus factor, we will pair up diagrams 
that differ by one encounter and one link to mutually can- 
cel before counting the surviving diagrams. The pairing 
is realized using the following recursive procedure: 

1. Find the outgoing leaf Om (attached to an encounter 
node) with maximal m. 

2. If Om is attached to a node of degree 4 whose op- 
posite edge ends in a leaf, untie the node. 

To untie a node we break it into two parts keeping the 
path connecting im to Om intact. This may separate a 
link directly connecting im to Om from the rest of the dia- 
gram. Such a link is removed from further consideration. 
We repeat steps 1 and 2 while it remains possible. For 
example, from Fig. [2fa) we separate the top encounter 
into two parts to give us Fig. [5Jb) and reduce m to 1. 

Once we can no longer perform step 2, we perform 
either of the following steps: 

3. If Om is attached to a node of degree 4 whose op- 
posite edge ends in another encounter node, shrink 
the edge and join the two nodes together. 

4. Otherwise, separate Om and its two neighbors from 
the encounter by inserting a new link. 

These two operations are inverses of each other and pro- 
vide the required difference in the number of encounters 
[see Figs. [21b) and (c)]. Finally we: 

5. reverse all the operations performed at step 2. 

We thus reconstruct a diagram paired to the original one, 
with a contribution of opposite sign. All diagrams that 
ever arrive at step 3 or 4 cancel with their partner, as, 
for example, the diagrams in Figs. Ufa) and (d). 

Factorizations of permutations. Diagrams that never 
reach step 3 or 4 can only involve encounter nodes of 
degree 4 and, following steps 1 and 2 repeatedly, must 



eventually end up as a set of independent links con- 
necting each ij to its Oj. For M„ we initially have tra- 
jectories along the boundary walk visiting the channels 
«i — ?• oi i2 . . . On — > ii which we can represent as the 
cyclic permutation (T„ = (12 . . .n). For systems without 
time reversal symmetry, when we arrive at step 2 for the 
first time we must have some Oj opposite o„. The oper- 
ation of untying is equivalent to multiplying (T„ by the 
transposition (j n) , breaking the boundary walk into two 
cycles, cr„(j n) — (1 . . . j)(j -I- 1 . . . n). Repeating steps 
1 and 2 we multiply repeatedly on the right by the pair 
of o channel labels at each step 2 until we arrive at the 
independent links whose boundary walk is the identity 
permutation. Reversing the untying of the nodes, we 
obtain a factorization of (t„ in terms of transpositions, 
which represents the original diagram. This process is 
illustrated in Fig. [31 Because we always chose the Om 
with the maximal m at each step 1, the resulting factor- 
ization {siti) . . .(sdtd) can be written so that tj > Sj 
and tk > tj for all k > j. Such a factorization is called 
'primitive' and its depth d is the number of nodes un- 
tied at step 2, each of which removes two links and one 
encounter. The number of encounters in the diagram is 
equal to d and the number of links n+2d. Diagrams with- 
out time reversal symmetry that survive the cancellation 
are therefore labeled by primitive factorizations and pro- 
vide the semiclassical contribution {—l)'^NiN2/N"^'^. 

Encounters in the lead. One complication is encounters 
that occur in the leads. For example, we could push 
the encounter in Fig. [ijb) to the left into the lead so 
that the incoming channels coincide ii — i2. We then 
lose two links, the encounter itself and one channel so 
that the new diagram still contributes at the same order 
but as NiN^/N^ instead of -N^N^/N^. Likewise we 
could move the encounter to the right until the outgoing 
channels coincide and the full contribution of diagrams 
related to Fig.[TJb) is the sum of these three possibilities. 

Whether encounters can be placed in the lead can be 
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(124){1342) 
={13)(2)(4) 




FIG. 4. (Color online) The top node in (a) may move into 
the outgoing lead which, by untying, we can represent as the 
permutation (1234)(124)=(1342) and the boundary walk in 
(b). Also placing the central node in (a) in the incoming lead 
means we further multiply on the left by (124) to obtain (c). 



seen directly from the graphical representation. If every 
alternate link of a node ends in an i leaf, then the en- 
counter can be placed in the incoming lead and similarly 
for the outgoing lead. Once encounters are in the leads, 
they can no longer be joined or separated as in steps 3 
and 4. Instead, we consider the encounters in the leads 
as being already untied. Then we perform the procedure 
above to identify canceling pairs. Each node in the out- 
going lead corresponds to multiplying (t„ on the right by 
the cycle of the labels of the o leaves (read off clockwise). 
Nodes placed in the incoming lead multiply the permu- 
tation on the left by the cycle of their i labels, illustrated 
in Fig. m For the resulting permutation all non-canceling 
diagrams are again labeled by the primitive factorizations 
of the permutation in question. 

Equivalence with RMT. The semiclassical result is 
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(1) 



where Pd{-) is the number of primitive factorizations of 
depth cZ, and c(-) the number of cycles. In the sum over 
all possible permutations Xi and Xo, their cycles represent 
incoming and outgoing channels which coincide. For ex- 
ample, the diagram in Fig. SJa) contributes to the fourth 
moment when the top and central nodes are placed in 
the leads, Xi = Xo — (124) (3). The resulting primitive 
factorization (13) leads to the contribution — TVf A/'l/iV^. 
For A/2, the only primitive factorizations are (12)'^ with 
odd and even d being factorizations of (12) and (1)(2) 
respectively, which are also the two possibilities for Xi 
and Xo- This gives M2 = N1N2/N - NfNi/{N^ - N). 

The random matrix result, on the other hand, 
can be written exactly as in Eq. ([T]) but repla cing 
Xld(-l)''pd(7r)/A^"+'' by coefficients V{Tr); see [s^ Efl 
41|. That V{t:) is the generating function for the number 
of primitive factorizations Pd{T^) was recently established 
[42| using an expression for ¥{11) in terms of characters 
of the symmetric group. Here we sketch a simple alterna- 
tive proof that is easy to generalize to systems with TRS. 
The Pdi^^) only depend on the cycle structure of tt so let 



ci , . . . , Cfc be the lengths of the cycles in the permutation 
TT = (12 . . .Ck) ■ ■ ■ {n — ci + 1 . . .n). First consider the 
case when the term on the right of a primitive factoriza- 
tion of TT is of the form [sd n). Without this term, it is a 
factorization (of depth d—\) of the permutation ■n{sd n). 
If Sd belongs to the rightmost cycle of tt it splits into two, 
of lengths q and r with g + r = ci , while if Sd belongs to 
cycle j the rightmost cycle joins with it to form a cycle 
of length ci + Cj. Finally, the last term (sdtd) can have 
td ^ n only if ci = 1 and the factorization does not have 
any occurrence of n in it. In this case it is also a factor- 
ization of the permutation on n — 1 elements with cycle 
lengths C2 , . . . , Cfc . In total we have 



PdiCl, . . . ,Ck) = (5ci,lPd(c2, 
+ '^CjPd-lici+Cj,...)-\- 



■■,Ck) 

^ Pd-i{q,r,C2,...,Ck) 



exactly mirroring the recursion relations of V [40| . 

Time reversal symmetry. With TRS the semiclassi- 
cal diagrams are more complicated , and the relevant 
combinatorial objects are permutations on 2n elements 
fi, . . . , 1, 1, . . . , n. The starting permutation with no en- 
counters in leads is encoded by (t„ = (n . . . 1)(1 . . . n). 
The cancellation procedure remains the same, but now it 
is also possible to untie nodes which have an ij leaf op- 
posite o,„, which we represent by multiplication on the 
right by (j m) and on the left by {fh j) [untying Oj and Om 
also multiplies by (mj) on the left in addition to {j m) on 
the right]. The non-canceling factorizations are then of 
the form pd(o'n) = [hsd) ■ ■ ■ {iiSi){siti) ■ ■ ■ {sdtd) with 
tj > Sj and tk > tj, tj, Sj for all k > j. 

Since the two leads for the transmission are separate, 
still only nodes with alternating i or o leaves can move 
into the leads and we obtain a result like Eq. ([T]) but 
involving the doubled permutations. Following similar 
reasoning to above, the Pd{T^) satisfy the same recursion 
relations as the coefficients V from the circular orthogo- 
nal ensemble so the moments are identical. 

Spin-orbit interaction. The semiclassical framework in- 
cludes spin-orbit interaction through an additional trace 
of a product of spin propagators along the trajectories 
[isf . The structure of the leading order diagrams makes 
this product identity and effectively leaves the leading 
order of M„ unchanged [S] . For each order higher in in- 
verse channel number, the chaotic spin-orbit interaction 
provides an additional factor of —1/2 ("for spin 1/2) com- 
pared to the contributions with TRS [44l|. The cancella- 
tion procedure then still holds, and the effect can be in- 
cluded by simply substituting A^i -> -2Ni, N2 -2N2 
and multiplying by —1/2. This is the same mapping as 
between the orthogonal and symplectic RMT ensembles 
[30I so the moments are again identical. 

Conclusions. Though we focused on the moments of 
the transmission eigenvalues, the combinatorial treat- 
ment here extends to non-linear statistics by simply 
changing the starting permutation ct„. Moreover, since 
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we have an exact concordance between the semiclassical 
and RMT [Hi^ moments of the transmission eigen- 
values, we obtain their probabihty distribution semiclas- 
sically. Indirectly, any RMT result derived from this 
distribution, for example the nonlinear statistics in 32 1 
and the moments of the conductance and shot noise 
3J, |45|, |46[, is now rooted in the chaotic dynamics in- 



side the cavity and the correlations between scattering 
trajectories. 

For energy dependent correlation functions, diagrams 
related to each other by the cancellation procedure no 
longer cancel exactly. For the related quantities of An- 
dreev billiards and Wigner delay times, the agreement 
between semiclassics and RMT remains limited to lead- 
ing [13, Ei] and several subleading orders p2l, [sl liij 
respectively. 

Finally, with the close connection between ballistic 
chaotic systems, RMT and systems with weak disorder, 
the combinatorial ideas here should have parallels in the 
diagrammatic perturbation theory of disorder. They also 
form the start point of including Ehrenfest time effects 
beyond the RMT regime. 
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